Anisotropic phase diagram of the frustrated spin dimer compound Ba 3 Mn 2 8 
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Heat capacity and magnetic torque measurements are used to probe the anisotropic temperature- 
field phase diagram of the frustrated spin dimer compound BasMi^Os in the field range from OT 
to 18T. For fields oriented along the c axis a single magnetically ordered phase is found in this field 
range, whereas for fields oriented along the a axis two distinct phases are observed. The present 
measurements reveal a surprising non-monotonic evolution of the phase diagram as the magnetic 
field is rotated in the [001]- [100] plane. The angle dependence of the critical field (H c i) that marks 
the closing of the spin gap can be quantitatively accounted for using a minimal spin Hamiltonian 
comprising superexchange between nearest and next nearest Mn ions, the Zeeman energy and single 
ion anisotropy. This Hamiltonian also predicts a non-monotonic evolution of the transition between 
the two ordered states as the field is rotated in the a-c plane. However, the observed effect is 
found to be significantly larger in magnitude, implying that either this minimal spin Hamiltonian is 
incomplete or that the magnetically ordered states have a slightly different structure than previously 
proposed. 

PACS numbers: 75.30.Kz, 75.40.-s, 75.30.-m, 75.45.+j 



I. INTRODUCTION 

Ba3Mri208 is a novel layered spin-dimer compound, 
comprising magnetic dimers of Mn 5+ ions arranged on 
triangular planes^. The Mn ions occupy equivalent sites 
with a distorted tetrahedral coordination, resulting in a 
quenched orbital angular momentum and a total spin 
S = l 2 . Antiferromagnetic intradimer exchange leads 
to a singlet ground state with excited triplet and quin- 
tuplet states. Weaker interdimer exchange leads to dis- 
persion of these excitations, as previously revealed by in- 
elastic neutron scattering 3 . The rhombohedral R3m crys- 
tal structure of BasM^Og comprises staggered hexago- 
nal planes (Fig. DJa)), leading to geometric frustration 
both within individual planes, and also between adja- 
cent planes. Magnetic fields can be used to close the spin 
gap, and the competition between interdimer interaction 
on this highly frustrated lattice, and the uniaxial single 
ion anisotropy associated with the S = 1 ions, leads to a 
very complex phase diagram with at least three distinct 
ordered states^*. 

Despite the complexity of the phase diagram of 
Ba3Mn208, we have previously shown that for fields ori- 
ented along the principal crystalline axes the high-field 
behavior of this material can be described by a remark- 
ably simple spin Hamiltonian, including terms repre- 
senting superexchange between nearest and next-nearest 
Mn ions within and between planes, the Zeeman en- 
ergy and uniaxial single ion anisotropy^. These terms 
have been determined through a combination of inelas- 
tic neutron scattering (INS) and electron paramagnetic 
resonance (EPR). INS studies revealed the dominant ex- 
change within a dimer as Jo = 1.642(3) meV; the nearest 
and next nearest out-of-plane exchanges, J\ = 0.118(2) 
meV and J4 = 0.037(2) meV; and the dominant in-plane 



exchanges J 2 - J 3 = 0.1136(7) meV (Fig. [Q) 3 . EPR mea- 
surements in the diluted compound Baa^i-^Mn^Os 
revealed a nearly isotropic g-tensor, with g aa = 1.96 
and g cc = 1.97, and an easy axis single ion anisotropy 
D = -0.024 meV— . Similar measurements in the pure 
BasMn2 08 compound revealed a zero field splitting of 
the triplet states characterized by D =- 0.032 meV— , the 
difference reflecting the additional effect of intradimer 
dipolar coupling in the undiluted compound. 

At low temperatures, the field dependence of the mag- 
netization of Ba3Mn208 reveals two regions of linearly 
increasing moment as first the S z = 1 triplet (between 
H cl = 8.7 T and H c2 = 26.5 T) and then S z = 2 quin- 
tuplet (between H c3 = 32.5 T and H cA = 47.9 T) states 
are polarized, separated by a plateau at half the satu- 
ration magnetization^. Here we concentrate solely on 
the first region of increasing magnetization, which can 
be accessed by moderate laboratory fields. In this field 
range the ground state can be approximated as a coher- 
ent mixture of singlet and S z =l triplet states. Previous 
measurements revealed a distinct anisotropy in the phase 
diagram depending on the field direction^. For fields ap- 
plied along the c axis, a single ordered state was found 
in this regime (fig. Gib) ). Analysis of the minimal spin 
Hamiltonian suggests that this state is an incommensu- 
rate XY antiferromagnet with wavevector shifted slightly 
away from the 120° structure favored for a 2D triangu- 
lar layer— NMR measurements appear to confirm this 
conclusion^, but to date the magnetic structure has not 
been explicitly solved. Alternatively, fields applied per- 
pendicular to the c axis revealed two ordered states in 
this field regime (Fig. [He)). Analysis of the minimal 
spin Hamiltonian for this field direction yields two dis- 
tinct incommensurate modulated states - an Ising phase 
stabilized closed to H c \ and iJ C 2, as well as a canted XY 
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FIG. 1: (Color online) (a) Schematic diagram illustrating 
the exchange interactions between Mn ions (red spheres) in 
BasMn208 (Ba and O ions are not shown). Values for the 
exchange constants are given in the main text (b, c) Phase di- 
agram of singlet-triplet ordered states for fields aligned along 
the c and a axes, respectively, taken from our previous work 4 . 
I, II and PM mark phase I, phase II and the paramagnetic 
phase as described in the main text. Solid and open symbols 
represent data points obtained via heat capacity and MCE, 
respectively. 

state in between. The modulation in these two phases 
is stabilized by the uniaxial anisotropy associated with 
the zero field splitting of the triplets. Similar to the sit- 
uation for fields oriented along the c direction, the exact 
form of the magnetic structure has not been explicitly 
solved, and it remains to be seen whether the minimal 
spin Hamiltonian that has been used to describe this sys- 
tem so far is complete. 

To this end, careful measurements of the phase di- 
agram of Ba3Mn208 for fields at intermediate angles 
between the a and c axes have the potential to deter- 
mine the presence or absence of additional terms. In 
this paper we present results of heat capacity and torque 
magnetometry measurements revealing how the two dis- 
tinct ordered states for fields perpendicular to the c axis 
evolve into a single phase for fields along the c axis. 
Through analysis of the previously established minimal 
spin Hamiltonian we can quantitatively account for the 
angular dependence of R c \ solely via consideration of the 
triplet dispersion. However, the same analysis, incorpo- 
rating the predicted magnetic structures, fails to quanti- 
tatively account for the angular dependence of the tran- 
sition between the two ordered states. We discuss the 
implications of this observation. 

II. EXPERIMENTAL METHODS 

Single crystals of BaaM^Os were grown from a NaOH 
flux according to our previously published procedure^. 




T(K) 

FIG. 2: (Color Online) Representative heat capacity data 
taken at 12T for for fields in the a-c plane. Labels indicate 
the angle between the field and the c axis. Successive data 
sets are offset vertically by 1.6 J/molK for clarity. 



Heat capacity (C p ) studies were performed on a Quantum 
Design physical properties measurement system (PPMS) 
using standard thermal relaxation-time calorimetry. 
These measurements were performed in fields up to 14T 
and temperatures down to 0.35 K. The sample was 
mounted on angled brackets made from oxygen-free high 
conductivity copper and the field was oriented in the 
[100]- [001] plane. 

Cantilever torque magnetometry experiments were 
performed at the National High Magnetic Field Labora- 
tory (NHMFL) in a superconducting magnet for fields up 
to 18T in a dilution refrigerator. A crystal was mounted 
on one face of a capacitance cantilever which was at- 
tached to a rigid plate rotatable about an axis parallel to 
the torque axis and perpendicular to the magnetic field. 
The magnetic field was aligned away from the principal 
crystalline axes yielding a finite torque. 



III. RESULTS 

Representative heat capacity measurements, taken at 
12 T for several angles, are shown in Figure [2j These data 
show a single peak for fields aligned along the c axis, a 
peak and a shoulder for fields 15° from the c axis and 
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FIG. 3: (Color online) (a) Phase diagram showing the tran- 
sitions between the paramagnetic state and phase II (Tjj), 
and between phase II and phase I (Tj), as a function of tem- 
perature and angle in the [100]- [001] plane for various fields, 
where indicates the angle between the field and the c axis, 
(b) Width of phase II, At = T CII — T C/ , as a function of an- 
gle for 11T, 12T, 13T and 14T (black circles, red up triangles, 
green down triangles, and blue squares, respectively.) 



two peaks for larger angles. Significantly, comparison of 
the data at 75° and 90° degrees shows that the 75° data 
has both a slightly higher critical temperature between 
the paramagnetic phase and phase II (T ClI ) and also a 
substantially lower critical temperature between phase II 
and phase I (T CJ ) than the 90° data. 

The phase diagram derived from the complete set of 
C p measurements, shown in Figure [3ja), reveals the evo- 
lution as a function of angle of the two distinct singlet- 
triplet ordered states for fields in the [100]- [001] plane. 
The data show a single transition for all fields for H\\c 
and two transitions for H\\a. The extent in temperature 
of phase II, At = T ClI — T CJ , is shown in Figure [3fb). 
At increases as a function of angle as the field is rotated 
away from the c axis, reaches a maximum at 75°, and 
decreases at the a axis (90°). For example, for a field of 
11T, A T is - 0.07 K larger at 75° than at 90°. 

Torque magnetometry measurements, taken at 25 mK, 
are shown in Figure |4] for three representative angles. 
The sample was inclined slightly so that the field did 
not exactly rotate within the a-c plane such that a finite 
torque was generated for all angles studied. Angles are 
quoted in terms of the angular position with respect to 
the closest approach to the c axis, but it is important 
to note that the field was never less than ~ 10° from 
the c axis. Consequently two phase transitions are ob- 
served for all angles studied. Critical fields, marking the 
transition between the paramagnetic phase and phase II 
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FIG. 4: Torque scaled by field, and its first two derivatives 
with respect to field, plotted versus field for representative 
angles at 5.5°, 55.0° and 88.0° with respect to the nearest 
approach to the c axis (see main text). The critical fields 
were estimated from peaks in the 2 nd derivative as shown by 
the dashed lines. 



(Hci), and between phase II and phase I (iJu_j), were 
determined from maxima and minima in the second field 
derivative of the torque divided by fields and are marked 
with dashed vertical lines. The sign of the peak in the 
second derivative changed with the evolution of angle, 
reflecting a change in anisotropy for the two different 
phases. This leads to minor discontinuities in the de- 
termination of the phase boundary (dashed lines in Fig. 

Ell- 
Similar to the phase diagram obtained from heat ca- 
pacity measurements (Fig. [3j), the phase diagram ob- 
tained from torque measurements at 25mK reveals a non- 
monotonic angle dependence (Fig. [5]). The maximum 
value of Hu-i occurs between 65-75° from the closest 
approach to the c axis. This is in agreement with the 
heat capacity data, for which the smallest T Cl occurs at 
75° from the c axis. Additionally, the field extent of phase 
II, Ah = Hi i -i — H c i, is largest at 75°, and decreases 
by ~0.05 T from the maximum at the highest angles. 



IV. DISCUSSION 

The previously established minimal spin Hamiltonian 
for arbitrarily oriented field direction in BaaM^Os is: 

u= J2 • S jv + D {SU cos - S£ sin 9) 2 

-HbH (9zzSfp + g X zS?p) , (1) 

i/xoi/3 



Numerical Calculation 



the critical field. The expression for H\\c is: 




FIG. 5: (Color online) Phase diagram at 25 mK determined 
from torque magnetometry measurements. Circles (triangles) 
mark transitions between the disordered phase and phase II 
(phase II and phase I). Open (closed) symbols signify that 
the transition was determined from a peak (trough) in the 
second derivative. Angles are measured relative to the closest 
approach to the c axis as described in the main text. Red 
dotted line shows the calculated transition between the para- 
magnetic and ordered phases as described in the main text. 



where g zz = g aa sm 2 + # cc cos 2 #, g xz = (g cc - 
g aa ) sintfcostf, g a p is the diagonal gyromagnetic tensor 
with components g CCl g aa = g^, and is the angle be- 
tween the applied field and the c-axis. The quantization z 
axis is set along the field direction. Here z, j designate the 
dimer coordinates, a,/? = {x,y,z}, v — {1,2} denote 
each of the two S=l spins in each dimer. The various ex- 
change constants are shown in Fig. QJa) and are denned 
as follows: the exchange within a dimer is Jo = Ji, 1,2,2; 
the dominant out-of-plane exchange is J\ = Ji,2,j,i fc> r 
z, j nearest neighbor dimer s between planes; the domi- 
nant in-plane exchanges between dimer s is J2 = Ji^j^ 
and J3 = Ji,/j,j,v for z, j in plane nearest neighbor dimers 
and fi ^ v\ and finally the second largest out-of-plane 
exchange is J4 = Ji,2,j,i fc> r hj next nearest neighbor 
dimers between planes. 

Using this spin Hamiltonian and the measured values 
of J0-J4 and D given in the introduction, we have pre- 
viously been able to quantitatively account for the ob- 
served critical fields and magnetization of the ordered 
states of Ba3Mn2 08 for fields oriented along the princi- 
pal axes 4 ^^. The calculation is based on a generalized 
spin-wave approach in which we only keep the singlet 
and the three triplet states of each dimer—. The criti- 
cal field H c i{0) corresponds to the value for which the 
energy of lowest energy triplet mode becomes equal to 
zero. The softening of this triplet mode signals the on- 
set of the magnetic instability towards an ordered state 
(phase II). For field directions along the principal axes, 
it is possible to obtain simple analytical expressions for 
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while for H JL c the expression is: 
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(3) 



where Jmin is the minimum of the interdimer exchange 
portion of the dispersion and is fully described in equa- 
tion (3) of reference [3[. The difference in the first two 
terms between these expressions stems from a change in 
the zero field splitting of an isolated dimer depending 
on quantization direction expressed in the reduced ba- 
sis of dimer states. For the quantization axis along the 
c axis, the zero field gap of an isolated dimer between 
the S z = 1 triplet and the singlet is Jo — D/3, while for 
the quantization axis along the a axis the equivalent zero 
field gap is Jo + D/6. The two additional terms of H c \ 
for fields perpendicular to the c axis arise from a sec- 
ond order process mixing singlets and triplets described 
in our previous work 4 . This state mixing causes the gap 
between the S z = 1 triplet and singlet states to close as 
\J H H c \ , as expected for an Ising-like QCP. Evaluat- 
ing these expressions using the values of the exchanges 
and single ion anisotropy described earlier yield values 
for the critical fields which are in good accord with the 
measured data. Numerical calculation of H c i for arbi- 
trary field orientations in the [001]- [100] plane yields the 
red dotted curve shown in Fig. [5] The calculation was 
performed using the values of D and the interdimer cou- 
plings given in the introduction while Jo was allowed to 
vary leading to a fit value of 1.567 meV. The calculated 
values agree well with the measured data up to the in- 
herent uncertainty associated with the misalignment of 
the sample in the torque measurements described above. 
The single-ion anisotropy term of H is: 



D (S z ) cos 2 (0) + (S x ) sin 2 (<9) 



(S Z S X + S X S Z ) cos (0) sin (0) 



(4) 



The last term of eq. (|4]) is zero for fields along the a 
and c axes, but adds a small contribution for interme- 
diate angles. In particular, this term grows linearly in 
small deviations of from 7r/2, 50 = ir/2 — 6, while the 
other two terms vary quadratically in 56. Thus, the last 
term of Eq. determines the shape of the boundary 
between phase I and phase II slightly away from = 7r/2, 
and could be responsible for the striking non-monotonic 
behavior observed for Hu-i in Fig. [5] as we describe in 
greater detail below. 

To understand the effect that the last term of eq. (|4j) 
has on the ground state, it is convenient to analyze the 
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effective low-energy Hamiltonian, H e ff, introduced in ref. 

The low-energy Hamiltonian results from projecting 
the original Hamiltonian T-L onto the low-energy subspace 
generated by the singlet and the S z = 1 triplet states of 
each dimer. This two- level Hilbert space is described by 
a local pseudo-spin ^ in each dimer, s^, such that s z 
if the dimer i is in the S z = 1 triplet state and s\ = - 
it is in the singlet state. In our earlier work we provided 
the expression of H e ff for = and = 7r/24 In par- 
ticular, for 6 = 7r/2, we showed that the second term of 
Eq. (j4j) generates an effective exchange anisotropy that 
is responsible for the emergence of phase II. According to 
our analysis, this exchange anisotropy favors an Ising-like 
phase in which the transverse spin components (perpen- 
dicular to the applied field) are aligned along the easy 
c-axis. In contrast, phase I is an elliptical spiral phase in 
which the transverse spin components of adjacent dimers 
rotate around the field axis. 

The effect the last term of Eq. (J2J) has on H e ff at in- 
termediate angles can be determined using second order 
degenerate perturbation theory. Such an analysis yields 
an effective Dzyaloshinskii-Moriya (DM) interaction be- 
tween dimers on adjacent bilayers connected by the J\ 
and J4 exchange constants: 




FIG. 6: (a) Schematic diagram illustrating the two dimers 
involved in processes of order J1D/J0 that leads to Eq. ([5]). 
The c axis is vertical, and the two dimers belong to adjacent 
planes. Ji represents either Ji or J4. (b) (i), (ii) Schematic 
diagram for two pseudo- spins on adjacent dimers connected 
by the Ji interaction for arbitrary field direction in the a-c 
plane. Red arrows show uniform moment, while green arrows 
show ordered moment. Field direction along z axis, Dz along 
y axis, and c axis vertical. 



Di-s^xsj-h ^2 D 4 -s^xsj, (5) 

with = Ay, Di = O(DJ t /J ) and I = 1, 4. The arrow 
indicates how the bonds ((i —> j)) are oriented (i always 
denotes the dimer in the lower bilayer). Microscopically, 
this process turns one singlet into an S z = 1 triplet or 
vice versa. 

This effective DM coupling between pseudo-spins 
results from two important symmetry considerations. 
First, the singlet and the triplet states of a given dimer 
have opposite parity under exchange of the two sites of 
the dimer: 1 <H> 2. Second, the J\ and J4 terms of H, are 
invariant under the inversion symmetry transformation 
around the center of the corresponding bonds: 1^2 
and i <B> j (see Fig. [6ja)). This implies that a DM 
term is allowed between pseudo-spins connected by the 
Ji and J4 exchanges. In contrast, the effective DM in- 
teraction cannot occur for pairs of dimers within a plane 
because the J2 and J3 terms of 1-L are invariant under 
the symmetry transformation i j, while the effective 
DM interaction changes sign under such transformation. 
These symmetries are fundamentally equivalent to the 
selection rules governing the DM vector in the effective 
lattice, where each dimer constitutes a single, unique site. 
In such a lattice, there is a center of inversion symmetry 
at the midpoint of the effective J2 and J3 exchange, pre- 
cluding a DM term by the normal selection rules, while 
no such inversion symmetry exists at the midpoints of 
the effective J\ and J4 exchanges. This is contrast to the 
symmetries of the real lattice, where there are inversion 
symmetries at the midpoint of the J\ and J4 exchanges 
but none at the midpoint of the J2 and J3 exchanges. 



This effective DM interaction between dimers on ad- 
jacent layers is frustrated in both ordered states at a 
mean field level (the mean value of Eq. (|5j) is zero for 
the semi-classical states associated with phases I and II). 
Therefore, the small contribution of the effective DM in- 
teraction to the ground state energy must be produced 
by quantum fluctuations. Because the DM vectors point 
along the y direction (perpendicular to the applied field 
and to the easy c axis) this contribution term will favor 
phase II, for which the pseudo-spins only have x and z 
components (Fig. [6fb)(ii)), as opposed to phase I for 
which the pseudo-spins have an additional third com- 
ponent along the y direction (Fig. EJb)(i))^. Thus this 
contribution strengthens phase II relative to phase I near 
= 7r/2 and leads to a small non-monotonic behavior of 
the Hi 1 -1 curve (see Fig. [5j). 

Although this simple analysis captures the qualitative 
non-monotonic behavior of the Hjj-j curve, it cannot 
account for the magnitude of the observed effect. The 
amplitude V{0) = Dsin^cos^ of the effective DM inter- 
action is of order 100 mK for c± 75°. Because the inter- 
action mixes the singlet and S z = 1 triplet dimer states, 
the mean value of the DM term is less than yJrnD(0) 
for any state with magnetization m (the mean density of 
S z = 1 triplets). Noting that m c± 0.06 at the 
upper bound on the DM term of c± 25 mK is the order of 
magnitude of the observed non-monotonic effect of 5-10 
mK in the Hu_j curve. Given that V(0) is much weaker 
than the dominant terms of it is clear that the ef- 
fective DM term can only explain the magnitude of the 
non-monotonic effect if it gives a first order contribution 
to the energy of Phase II. However, as established above 
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for the proposed ordered states, the effective DM interac- 
tion contributes via a second order correction and must 
therefore be considerably smaller. This leaves us with 
two possibilities: a) The magnetic structure of Phase II 
is different from the simple Ising phase proposed in Ref. 
0] in such a way that the mean value of the effective 
DM term is non-zero, or b) The non-monotonic effect is 
caused by a term that has not been included in At 
present it is impossible to distinguish between these pos- 
sibilities, but ongoing efforts to experimentally determine 
the magnetic structure have the potential to directly ad- 
dress option (a), while EPR experiments should, at least 
in principle, be able to determine the energy scale of ad- 
ditional interactions not considered in the minimal spin 
Hamiltonian (eq. (pQ)). 

V. CONCLUSION 

In summary, via heat capacity and torque magnetome- 
try measurements we have established the angular depen- 
dence of the phase boundary for singlet-triplet ordered 
states of Ba3Mn2 08- The data reveal a striking non- 
monotonicity of the phase boundary as the field is rotated 
between the principal axes. The angle-dependence of H c \ 



can be quantitatively understood in terms of the origi- 
nal minimal spin Hamiltonian that we had proposed for 
this material. This quantity does not depend on details 
of the magnetically ordered states but only on the min- 
imum of the triplet dispersion. However, the observed 
non-monotonicity in Hi-n is at least an order of mag- 
nitude larger than anticipated based on this model and 
assuming the magnetic structures previously proposed. 
This indicates that a complete theoretical description of 
this material requires either subtle changes in the pro- 
posed magnetic ordered structures or an additional low- 
energy term in the Hamiltonian. 
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